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T^his  is  a sequel  to  a previous  paper  where  rational  approxi- 
mations for  the  confluent  hypergeometric  function  zaU(a;c;z) 
'were  treated.  Here  we  take  up  rational  approximations  for 
^,F^(a;c;-z).  The  confluent  functions  are  very  important  in  the  • 
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applications  since  they  include  as  special  cases  the  incomplete 
gamma  function^ (special  cases  of  which  are  exponential,  sine  and 
/Cosine  integrals,  Fresnel  integrals  and  the  error  function), 
^Bessel  functions,  parabolic  cylinder  functions  and  Coulomb  wave 
functions . -^.The  subject  of  rational  approximations  for  a wide 


class  of  functions'  including  those  named  above  were  examined  in 


some  detail  in  my  volumes  on  the  special  functions . ^ In  the 
special  case  where  a is  unity,  the  confluent  function  becomes  an 
incomplete  gamma  function.  In  this  event,  complete  a priori  error 
analyses  for  the  main  diagonal  Fade  approximations  and  much  more 
were  presented.  For  general  parameters,  the  rational  approxima- 
trons  treated  were  not  of  the  Pade  class.  It  was  shown  that  the 
rational  approximations  converge,  but  a complete  a priori  analysis 
was  not  available.  One  of  the  purposes  of  this  report  is  to  cor- 
rect this  deficiency.^ Further , FORTRAN  programs  are  provided  to 


evaluate  the  Pade  and  non-Pade  rational  approximations  by  using 
the  appropriate  recursion  formulas  to  generate  the  numerator  and 
denominator  polynomials  as  a number,  and  to  also  evaluate  the 
coefficients  which  define  these  polynomials.  The  programming  of 
the  routines  was  done  for  use  by  the  IBM  370/168  operating  under 
OS/VS  Release  1.7  on  the  FORTRAN  IV  H-Extended  Compiler,  Release 
2.1.  All  computer  programs. are  written  for  quadruple  precision 
and  real  arithmetic.  By  making  a few  simple  changes,  one  can  have 
double  or  single  precision..  Further,  it  is  easy  to  get  complex 
arithmetic  along  with  any  of  the  precisions  noted  above. 
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Summary 

This  is  a sequel  to  a previous  paper  where  rational  approxi- 
nations  for  the  confluent  hypergeometric  function  zaU(a;c;z) 
were  treated.  Here  we  take  up  rational  approximations  for 
jF^(a;c;-z).  The  confluent  functions  are  very  important  in  the 
applications  since  they  include  as  special  cases  the  incomplete 
gamma  function  (special  cases  of  which  are  exponential,  sine  and 
cosine  integrals,  Fresnel  integrals  and  the  error  function), 

Bessel  functions,  parabolic  cylinder  functions  and  Coulomb  wave 
functions.  The  subject  of  rational  approximations  for  a wide 
class  of  functions  including  those  named  above  were  examined  in 
some  detail  in  my  volumes  on  the  special  functions.  In  the 
special  case  where  a is  unity,  the  confluent  function  becomes 
an  incomplete  gamma  function.  In  this  event,  complete  a priori 
error  analyses  for  the  main  diagonal  Pad^  approximations  and 
much  more  were  presented.  For  general  parameters,  the  rational 
approximations  treated  were  not  of  the  Pad£  class.  It  was  shown 
that  the  rational  approximations  converge,  but  a complete  a priori 
analysis  was  not  available.  One  of  the  purposes  of  this  report 
is  to  correct  this  deficiency.  Further,  FORTRAN  programs  are 
provided  to  evaluate  the  Pade  and  non-Pad£  rational  approximations 
by  using  the  appropriate  recursion  formulas  to  generate  the 
numerator  and  denominator  polynomials  as  a number,  and  to  also 
evaluate  the  coefficients  which  define  these  polynomials.  The 
programming  of  the  routines  was  done  for  use  by  the  IBM  370/168 
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operating  under  OS/VS  Release  1.7  on  the  FORTRAN  IV  H-Extended 
Compiler,  Release  2.1.  All  computer  programs  are  written  for  - 
quadruple  precision  and  real  arithmetic.  By  making  a few  simple 
changes,  one  can  have  double  or  single  precision.  Further,  it 
is  easy  to  get  complex  arithmetic  along  with  any  of  the  precisions 
noted  above. 
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1.  Introduction 

In  a previous  paper,  Luke  (1),  we  considered  rational  approxi- 
mations for  the  confluent  hypergeometric  function  z*U(a;c;z). 

This  is  a sequel  to  the  paper  just  noted  in  which  we  treat  rational 
approximations  for  the  confluent  function 

CD 
(2) 

The  connection  between  the  confluent  functions  is  given  by 

U(a;c;z)  - ^(ajcjz)  ♦ ^(b^-cjz) , b - l*a-c. 

(3) 

The  confluent  functions  are  very  important  in  the  applica- 
tions since  they  include  as  special  cases  the  incomplete  gamma 
functions  (special  cases  of  these  are  exponential,  sine  and  co- 
sine integrals,  Fresnil  integrals  and  error  functions),  Bessel 
functions,  parabolic  cylinder  functions  and  Coulomb  wave  functions. 
For  more  details  on  all  these  transcendents,  see  Luke  (2,3). 


Consider  a function  E(z), 

E(z)  - En(z)  ♦ Rn(z), 

(4) 

Bn(z)  " AnCO/Bn(z) 

(S) 

whete  Ar(z)  and  Bn(z)  are  polynomials  in  z of  degree  n. 


B(z)  - 1?1(a;c;-z), 


k k 


- (»)k(-r* 

B(s)  " 2 — fcT ""1 1 
k-0  (CJk** 


— 


" 


■n  !■>. 


mam  - hhi 
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and  (z)  is  the  remainder.  For  each  E(z),  with  parameters 
specified,  two  FORTRAN  programs  along  with  illustrative  examples 
are  given.  In  the  first,  a positive  integer  N and  z are 
also  specified.  The  machine  is  given  forms  to  evaluate  sufficient 
initial  values  of  Ar(z)  and  Bn(z),  and  computes  subsequent 
values  of  An(z)  and  Bn(z)  by  means  of  a recursion  formula  and 
evaluates  En(z)  for  n - 0,1, ...,N.  The  machine  prints  An(z), 
Bn(z),  En(z),  En+1(z)  - En(z)  (that  is,  first  differences)  and 
' En  (z) , all  for  n • 0,1,..., N.  The  last  two  quantities 
can  be  viewed  as  a measure  of  the  error  with  the  latter  preferred. 
This  program  is  especially  valuable  both  to  get  the  desired 
approximants  En(z) , and  in  the  absence  of  applying  the  asymptotic 
forms  of  the  remainder,  to  also  get  an  appraisal  of  the  value  of 
n needed  to  achieve  a given  level  of  accuracy.  Once  n is  known, 
it  might  be  more  convenient  and  economical  to  have  the  coefficients 
which  define  the  polynomials  in  An(z)  and  Bn(z).  This  is  fur- 
nished by  the  second  program.  Thus  once  z is  specified,  A^C*)* 
Bn(z)  and  En(z)  are  readily  found.  The  polynomials  are  espe- 
cially valuable  when  one  desires  to  use  the  approximation  to  sim- 
plify analytical  formulas  and  to  make  further' approximations  such 
as  the  evaluation  of  integrals  and  transforms  involving  E(z). 

The  technique  for  evaluating  the  coefficients  in  the  numer- 
ator and  denominator  polynomials  follows.  It  is  convenient  to 
treat  more  general  forms  than  is  required  in  the  present  analysis. 
Let 


B„U) 


Ln2  q+f+3Fp+g+l 


/ - n , n+  X , - a , c £ , 1 

\p+l,  <»p+l-a,  dg 


-1/z 


■ 


(6) 


s 


V*>  “ V 


n(n+X)  (pQ-l)cf: 

— (HT^r 


] 


a 


n-a  (a-n)k(n+X+a)k(a  )k(c£+a). 

* ^i«»)kvi>kV‘7?r 


* q+f+3Fp+g+l 


(-n+a+k,n+x+a+k,p^+k,C£+a+k,l  \ ^ 

0+1+a+k,  a +l+k,  d +a+k  ) 


or 


An(D 


n(n+X) (p  -l)cf 

-TFTjrr — 

p i 


Lnz 


n-a 


n-a  (a-n)k(n+X+a)k(pQ)k(cf+a)k(-)  z 

\io lH‘.)k|.p3i]k(V*)i 


k -k 


-n+a+k,n+X+a+k*pq+k,c£+a*k,Op 
p+q+f+2  p+q+g+1  | p+i+a+kt  ap+l+k,  dg+a+k,  p^ 


1 • (8) 


where 


Ln  " 


(p+l)_ (o  ♦i-a)  (d  ) 

(nA)°n^Pq-,^tcBf)nn  • • ’ • « • ’ I-  * ■ ~**' 


(9) 


It  is  easily  proved  that 

/- n-p,  -n-o  +a,  1-n-d 

Bn^*^  " p+g+lFq+f+l  ( . , . 

n p x * h y l-X-ln.l-n-Pq+a.l-n-c^ 


(-)T*  . 


• Z u.z  , 
k-0  * 


uQ  - 1,  r * p+q+f+g. 


(10) 
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Then  the  u^'s  are  readily  evaluated  by  use  of  the  recurrence 


formula 


(-)r(-n-  p+k)  (-n-a^+a+k)  (l-n-de+k)uk 
Uk+1  " (l-\-2n+k)  (1-n-p  +k+k)  U'-n-c*+fc)  (k+1)  * 


Further,  we  can  write 


A (z)  « 2 v,z  , v,  - 2 t u.  _ , 

n k-0  k k m-0  m k‘m 


t - lfiy* 


where  the  uk's  are  defined  by  (10)  and  are  easily  evaluated  by 
(11).  Since 

(a  +m) 

*■♦1  " ' (pq+m) (m+1)  tm»  to  = 1 » 

we  see  that  computation  of  vk  is  direct. 

In  Luke  (2,3),  a is  either  0 or  1,  but  in  our  present  work 
a ■ 0.  This  a has  nothing  to  do  with  the  a in  ^F^(a;c;-z). 

We  now  show  how  to  get  the  forms  in  Sections  2 and  3.  Put  X ■ p+1, 
p « q * 1.  To  get  the  Section  2 polynomials,  let  P“f»g-0, 

- a and  p^  * c.  To  get  the  Section  3 polynomials,  let 
P ■ c-1,  take  f ■ g ■ 1 and  put  ■ 1,  p^^  * c,  c^  - 2 and 
di  " 1* 

As  in  our  previous  paper,  the  programs  are  written  for  certain 
IBM  equipment  in  quadruple  precision  for  real  arithmetic.  However, 
with  a few  minor  changes,  we  can  have  single  or  double  precision. 
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and  for  any  precision  we  can  also  have  complex  arithmetic.  All 
of  this  is  well  detailed  in  Luke  (1,  Section  4)  and  will  not  be 
repeated  here. 

2.  Rational  Approximations  for  ^jCajcj-z) 

Let 

E(z)  - ^(ajcj-z)  - e'z  ^(c-ajcjz) . (1) 

We  suppose  that  neither  a nor  c-a  is  a negative  integer  or 
zero,  for  otherwise  from  (1),  E(z)  is  a polynomial  except  for 
the  possible  presence  of  e"z.  We  also  suppose  that  c t a.  This 
is  for  convenience  only,  because  when  c - a,  E(z)  - 0FQ(-z)  * e*z, 
and  in  this  event  one  should  use  the  rational  approximations  of 
Section  3 with  c ■ 1 as  they  are  simpler  and  more  accurate.  The 
Section  3 approximations  with  c - 1 can  be  found  from  those  of 
this  section  by  putting  c * a+1  and  letting  a -♦  «.  For  general 
parameters,  the  rational  approximations  which  follow  are  not  of 
the  Pad£  class.  We  write 


E(z)  = {An(z)/Bn(z)>  + Rn(z)  , 


(2) 


Bn(z)  * LnzI1  jFjC-n.n+l.cja+lj-l/z) 


'-n+k,n+l+k,c+k,l 


n (-n).(n+l)k(a), 

A„(*)  " Kz  z -7--  ,F7 

n n k°0  (a+1).  (k'.)z  4 1 \l+k,  a+l+k 


(3) 


1/z  \» 


C4) 
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(a+l)n 

Ln " (nV1Vc>n  ' (5) 

Here  Rn(z)  is  the  remainder  which  we  discuss  later. 

For  the  polynomials  Bn(z)  and  A (z) , we  have 

B„(z)  - 1.  Bl(z)  - 1 ♦ isjin.  , B2(z)  - 1 ♦ . 

fa+3)z  . (a+2) 2z  (a+l)3z3 


b3c^  * 1 + + i(j(d4)2  + -izo(o3  * 

A0Cz)  - 1,  Al(z)  - Bl(z)  - , A2(z)  - B2(z)  - *§  [l  ♦ £f$j] 


(a),z< 


+ it hit  ' a3(z)  " Vz>  - H 


. fa+3)z  _ (-a+2^2z2 

rcc^iy  nr(c+T32 


- (a)2Z  Tl  - (a+3)zl  (a)3z3 

ttctt  l1  ■ st^tt  • 


(6) 


Both  An(z)  and  Bn(z)  satisfy  the  same  recurrence  formula 

Bn(z)  - (1  ♦ F1z)Bn.1Cz)  ♦ (E  + F2z3zBn.2(z)  ♦ F3z3Bn_3(z) 


n " 3»  F1  * IT2n?jj(nJc-iJ  » F2 

(n+a-2)-(n-a-2) 

Cm—.  • 


8 (2n-3; (2n-5) (n+c-3). 


, E 


(n+a) (n+a-1) 

T (In  -1 ) (2n  - 5y'(h+c  - 2 ) , » 


Ln+a-jQ  In  - c-lj  . ^ ? ^ 

r-TjT 


-3) in+c- 


The  recurrence  formula  is  stable  in  the  forward  direction.  We 
write  the  error  in  the  form 


Vz>  • S„(z)/Bn*(z),  Ba*(z)  ■ *'V\(z) 


(«) 
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and  first  consider  the  structure  of  and  asymptotic  forms  for  Sn(z). 
A closed  form  representation  for  SR(z)  can  be  derived  after  the 
manner  of  discussion  given  by  Fields  (4-6).  We  have 

Sn(z)  - F(z)Mn(z)  + H(z)Gn(z),  (9) 

where 


Asymptotic  forms  for  f'n(z)  and  Gn(z)  for  n large  with  a,  c 
and  z fixed  follow  from  the  work  of  Luke  (2,3).  Thus 


(c,  c-a 
c-n  ,n+c+l 

[1.+  0 (n*2)  ] . 


0(n"2r"2) 


(13) 

(14) 


Gn(z) 


(-)n+1n!  r (n-'-l-a)  zn+1 

r(n+i-c5(Zn7lT! 


exp 


[1  vz 4 


0(n 


-2 


v 


(n+l-a) 

8 (2n+3) (n+3-c) (n+2-c)2 


n! r (n+l-a)  _ n!nc‘a‘1 

f (n+2-c) (2n+l) : TTnTTJT- 


[n'1  + n(a-3c+6)  + (a-5c+ac+7) ] . , 
[1  + Ofr'1)].  (15) 


If  c is  not  a positive  integer,  then 
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F C2) r (n+l-c)  cr 


c,  c-a 


-2r-2 


Vz>  " 2F2  I n+1+c  2 + °(n  ' * <16> 

\ c-n,n+l+c 

If  a is  not  a positive  integer,  but  c is  a positive  integer, 
then 


V»>  • «<»>  ,F|  r . * 0(n'2r-2)  , 

c-n ,n+l+c 


-2r-2. 


where  E(z)  is  given  by  (1).  If  both  c and  a are  positive 
integers,  c > a,  then 


Sn“> 


az  r (c-l)e  r,a-l 


rfc-a] 


Fx  (l-a;2-c;z)Gn(z). 


With  a,  c and  z fixed,  it  can  be  shown  from  the  work  of  Luke 
(2,3)  that 

VCO  ■ 7— 77— s~ 7 [«p  a - «2  * OC"'2)]  , 

(a+l)_z  n!  L 


„2  (n+c-1' 


8(2n-l) (n+c-1) (n+c-2) 


[n  +n{ (3-2c)(2-a-c)+2(c-a-l) )+ac+2-2c] 


If  the  numbers  c and  a are  arbitrary  except  as  previously 
noted  with  the  further  provision  that  if  these  numbers  are  positive 
integers,  they  are  not  so  simutaneously , then  the  forms  for  the 
error  readily  follow  from  (8),  (9)  and  (16)  or  (17)  as  appropriate. 
In  these  situations,  we  have 


M*)  - 


Wr (n+l-c) r (n+a+l)n! z 
F (n+l+c j r fn+c) (2n) ! 
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where  W is  free  of  n except  that  it  might  contain  the  factor 


(-)n-  Clearly 


and 


|Wz^RnCz)|  = 


lim  R (z)  = 0 
_ n v 
n-*— 


(n+l-c) (n+a+1) z 
2 (n+c)"(n+c+l)  (2n+l)' 


(21) 


[1  + 0(n_1)  ] . (22) 


If  both  c and  a are  positive  integers  with  c > a,  then 
from  (8),  (9),  (15)  and  (18)  we  have 

,n_2n+2-c_-zrW_AT,,_  , -v  <•_  i •>  2 


1-a 


v*>  ■ i^r1  u.e 


[i  * <u*v> 2 2 * °("‘2 


)].  (23) 


Again 


lim  R (z)  - 0 , 
n-*— 


(24) 


and 


Rn+l(z)/Rn(z)  " ' yan*l}(2K?j  fn^c)(n+2-c)  [1  + 0(n  1)] 


(25) 

It  is  of  interest  to  compare  the  error  R^(z)  where  a * 1 
with  the  corresponding  error,  call  it  Rr  p(z),  for  the  Pade 
approximations  in  the  next  chapter.  If  c is  not  a positive 
integer, 

R_  _(z)  r ^n+1  n z/2 , . .3/2  . 

-Hrr  " U ■ [1  + °Cn  * C26) 

RntzJ  r(c)T(c-l) (sin  trc)n!2Zn  Zc  1 


whence  the  Pade  approximation  is  superior.  Now  consider  (23). 
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If  z is  small, 

Rn(z)  - 0(z2n+2_c).  (27) 

This  would  be  the  situation  for  a Pade  approximation  if  the  numer- 
ator and  denominator  polynomials  were  of  degree  n+l-c  and  n 
respectively.  Thus  our  rational  approximations  under  the  condi- 
tions leading  to  (23),  though  not  of  the  Pade  class,  are  very 
much  akin  to  this  class.  Indeed,  when  a ■ 1 and  c is  a posi- 
tive integer,  c > 1 , we  find  that 

- - (z/4n)c  1ll  + 0(n  X)]  , (28) 

and  again  the  Pade  approximation  is  superior. 

Numerical  Examples 

Let  a ■ 2/3,  c ■ 4/3,  z = 3/4  and  n ■ 5.  Then  from  (8), 

(9)  and  (16)  without  order  term? , 

Bn(z)  - 1.47509,  Bn*(z)  - 1.04541(5),  Rn(z)  - -0. 23921(-7) . 

The  true  values  of  Bp(z)  and  Rn(z)  are  1.47779  and  -0. 23914 (-7) , 
respectively. 

Suppose  a ■ 1,  c * 2 and  z * 5/4.  If  n * 5,  then  from 
(23)  without  order  terms,  Rn(z)  ■ -0.2762(-9)  which  agrees  with 
the  true  value. 
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3.  Rational  Approximations  for  jFj(l;c;-z) 


Let 


E(z)  - 1F1(l;c;-z) 


(1) 


which  is  a form  of  the  incomplete  gamma  function.  See  Luke  (2,3). 
The  rational  approximations  in  this  section  lie  on  the  main  diagonal 
of  the  Pade  table.  If  c * 1,  we  get  approximations  for  the 
exponential  function,  see  later  comments. 

We  write 


E(z)  - (An(z)/Bn(z)>  + Rn(z)  , 


n n (-n).(n+c), 

V*>  ' V kJ0  — TcT^l — sFi 


-n+k,n+c+k,l 

1+k 


,1 

z 


. r(n+c) 
-n  fffn^cl 


(2) 

(3) 

(4) 

(5) 


where  Rn(z)  is  the  remainder  which  is  discussed  later. 

For  the  polynomials  An(z)  and  Bn(z),  we  have 

B0(z)  - 1,  B1(z)  - l + , b2(z)  - l + ♦ T^Tr-  , 

b3(2)  " 1 + Hr  + Tc^tyj  + TcUyi  •* 

Aq(z)  - 1,  Ax(z)  - Bx(z)  - | , A2(z)  = B2(z)  - |(l  + + j|y- 


r 
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a3(z)  - Bj(z)  - f(i  ♦ §fy  ♦ x^rfyr  ♦ xfr:!1  + Hr)-  xItT  * (6) 


Both  A^z)  and  Bn(z)  satisfy  the  same  recursion  formula 
Bn(z)  - Cl  ♦ Fiz)Bn.1(z)  + F22ZBn.2(i),  n > 2 , 

„ fc-1)  „ _ fn-1) fn+c-2) 

F1  (Zn+c-15 lin+c-S)  • *2  (2..e-l)(lm*e-S)‘C2..e-«)  ' (?) 


The  recursion  fornula  is  stable  in  the  forward  direction. 
For  the  remainder,  we  have 


v*>  * 


f-)n*^nr(c)n!r(n»c)z^n+1{exp  [*z+z(z+4c-4)/4(2n+c)])  [l*0(n  3)), 


24n+2c-2(2n+c)  [p(n+c/2)r (n+ (c+l)/2) ] 


or 


Vz> 


(-)n+1Trr(c)nx'cz^nr x{exp  [*z+z(z+4c-4)/4(2n+c)  3)  [l+0(n  x)3 


l-c_2n+l 


2*4n+2c-i  (n!)2 


(8) 


-1, 


It  follows  that  for  z and  c fixed, 

lim  R fz)  * 0. 
n 

n-*« 


(9) 

(10) 


To  facilitate  computation  of  a priori  error  estimates,  we 


have 


Rn+1^ 


”P  [2(2^c)(2n*c5Tr]  C1*0(n'S” 


(2n+c) (2n+c+l)  (2n+c+2) 
.2 


z"  j [1  ♦ 0(n"2)]  , 


4(2n+c+l) 


(ID 


which  is  a measure  of  the  rate  of  convergence.  With  Rft(z,c)  in 


A 


IS 


place  of  Rn(z) , 

-\f0  £^THT1  (tH+c)  (2n*c^l)  ) exP  [Tr7nicH2n*c*h)] 

(12) 

The  latter  shows  that  for  a given  z,  c and  n,  the  error  changes 
but  slightly  for  snail  values  of  h,  that  is,  for  snail  changes  in 
c. 

Next  we  consider  sone  results  for  the  exponential  function  - 
the  case  c ■ 1.  We  have 


•"*  - {Gn(-*)/Gn(z)>  ♦ Sn(z)  , 


(13) 


where 


(14) 


Gn(z)  - z11  2P0(-n,n+l;-l/z)  - -^1  ^(-n; -2n;z)  , (15) 


«nd  Sn(z)  " ^n^z^  with  c ■ 1.  It  is  convenient  to  write 

Gn(z)  “ Mn(z2)  ♦ *Nn(*2)* 


(16) 


Then  by  conputing  ^(z2)  and  Nn(z2),  evaluation  of  the  nain 

diagonal  Padl  approxination  of  order  n only  necessitates  the 

evaluation  of  essentially  (n+1)  terns.*  The  polynonials  Gn(z), 
2 2 

Mn(z  ) and  Nn(z  ) satisfy  the  same  recurrence  fornula. 


Gn+l(z)  " 2(2n*l)Gn(z)  ♦ z2Gn.1(z)  , (17) 
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GQ(z)  - 1,  GjCz)  - z+2 , G2(z)  - z2  + 6z  ♦ 12  , 
G3(z)  - z3  ♦ 12  z2  ♦ 60  z + 120, 

G4(z)  ■ z4  + 20  z3  + 180  z2  + 840  z + 1680. 

We  also  have  the  explicit* representations , 


%n,  % - %n 

2 \ 

z2/4 

» 

n,  h - n, 

/ 

Ih  - ta,  1 ■ 

■ %n 

2 

*V 4 

\1  - n,  % - 

n,  3/2 

N„(z2)-^?^  2F3P,n''i] 

where  [y]  stands  for  the  largest  integer  in  y. 

As  previously  noted,  forms  for  the  error  sn(z)  follow 

(8)  - (11)  with  c ■ 1.  In  closed  form, 

(-)n+1*e-zI  (z/2) 



and  we  have  the  asymptotic  representation 

Sn(z)  - (-)n+1e‘z[exp(2vc  ♦ 2U/v)][l  ♦ 0(v‘3)] 

uniformly  in  z,  z * 0,  |arg  z|  < tr/2,  where 

v - n ♦ >j,  z • 2vx , c - u"1  ♦ , u - (1  ♦ x2)'1*. 


(18) 

(19) 

(20) 

from 

(21) 

(22) 


U - (3u  - Su3)/24. 


(23) 
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In  particular,  for  z large,  we  have 

Sn(z)  - (On+1  exp Qj  (l  - * 0(x‘3))]  exp  [^]  [l+0(v'3)]. 

(24) 

In  illustration,  let  n - 4,  z ■ 9 whence  v * 9/2.  Neglecting 
order  terms,  from  (22)  we  get  S4(z)  ■ -0.01474  whereas  the  true 
value  is  -0.01503.  Uniform  asymptotic  representations  for  the 
first  and  second  subdiagonal  Pade  approximations  for  e'z  are 
also  given  in  (7).  For  further  comments  on  Pade  approximations 
for  (1)  and  other  remarks  on  the  exponential  function,  see  (2,3,7, 
8,9).  In  connection  with  reference  (7)  Dr.  M.G.  de  Bruin, 
Universiteit  van  Amsterdam,  Instituut  voor  PropedeUtische  Wiskunde, 
Roetersstraat  15,  Amsterdam,  Netherlands  has  kindly  informed  me 
that  in  quoting  the  Tesults  of  H.  van  Rossum,  I overlooked  the 
restriction  u < v+1.  Consequently  the  general  results  of  refer- 
ence (7)  are  not  valid  for  the  lower  part  of  the  Pade  table  unless 
c ■ 0.  The  case  c * 0 is  for  the  exponential  function.  Also 
the  concept  of  'normal*  employed  by  H.  van  Rossum  is  different 
fTom  that  usually  employed  for  the  Pad£  table. 

Numerical  Examples 

Let  c ■ 2.  Then 
1F1(ls2;-x)  - (1  - e'z)/z. 

-We  take  z * Tf " V.  (n)  is  the  right  hand  side  of  (8)  with 


r 
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-3 


0(n  ) neglected,  then 


.n+1 


Vi(n)  - LT  x exp[-(16n,r7lZ32(nnj]  # 
1 2vii^Z  r/.,  tZ 


C(3/2)n]' 


Values  of  V^n)  and  the  true  values  of  V(n)  for  n • 0(1)3 
aTe  recorded  in  the  following  table. 


(-)n+1V1 (n) 


(-)n+1V, (n) 


0 

1 

2 

3 


0.201 
0. 121  (-2) 
0. 2 8 9 ( - 5 ) 
0 . 360 (-8) 


0.213 
0. 122 (~2) 
0.290(-5) 
0. 361  (-8) 


Thus  the  values  of  (n)  are  remarkably  accurate  even  for  small 


values  of  n.  If  V2(n)  is  the  right  hand  side  of  (9)  with 


-1 


0(n  ) neglected,  then 


V, (n)  - t~^n  1 ex: 


Pn[-jl6n»7)/32(n»1)] 


Values  of  V2(n)  and  the  true  values  of  V(n)  for  n ■ 3,4,5 
are  recorded  in  the  following  table. 


t-)n4lv2(») 


3 

4 

5 


0.4Sl(-8) 
0. 326 (-11) 
0.161 (-14) 


0.361(-8) 
0.274 (-11) 
0. 140(-14) 


The  above  is  the  sample  problem  treated  in  the  programs  which  follow. 
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THIS  MAINLINE  PROGRAM  ILLUSTRATES  THE  SYNTAX  OF  THE  SUB- 
ROUTINE ' RlFl ' FOR  GENERATING  VALUES  OF  THE  NUMERATOR  AND 
DENOMINATOR  POLYNOMIALS  IN  THE  RATIONAL  APPROXIMATION  OF 
IF1(  AP  ; CP  ; -Z  ) . 


IMPLICIT IUAL*16 (A-H , O-Z) 

DIMENSION  A(26) ,B(26) ,R(26) ,D(26) ,E(26) 
DATA  ZERO/0. QO/ 

R£AD(5,1, END-999)  N,M 
READ(5,2).AP,CP. 

Nl-N+1 
D(N1)-ZER0 
E(NI)-ZERO 
DO  100  I-1,M 
READ(5.2)  Z 


CALL  RlFl  (AP , CP , Z , A,  B , N) 


IN  THE  ABOVE  : 

AP  IS  THE  NUMERATOR  PARAMETER  OF  THE  1F1 

CP  IS  THE  DENOMINATOR  PARAMETER  OF  THE  1F1 

Z IS  THE  VALUE  OF  THE  ARGUMENT 

A AND  B WILL  CONTAIN  THE  VALUES  OF  THE  NUMERATOR  AND  DENOMINATOR 
POLYNOMIALS,  RESPECTIVELY,  FOR  ALL  DEGREES  FROM  0 TO 
N INCLUSIVE 

N IS  THE  MAXIMUM  DEGREE  FOR  WHICH  VALUES  OF  THE  POLYNOMIALS 

ARE  TO  BE  CALCULATED 

NOTE  : VALUES  OF  THE  K-TH  DEGREE  POLYNOMIALS  WILL  BE  PLACED  IN 
A(K+1)  AND  B(K+1)  RESPECTIVELY. 


R(N1)-A(N1)/B(N1) 

DO  50  J-l.N 
Jl-Nl-J 

R(J1)-A(J1)/B(J1) 

D(J1)-R(J1+1)-R(J1) 

50  E(J1)-R(N1)-R(J1) 

WRITE (6, 3)  N,AP,CP,Z 
DO  60  J-l.Nl 
Jl-J-1 

60  WRITE(6,4)  Jl,A(J),B(J) 

WRITE (6, 5) 

DO  70  J-1,N1 
• Jl-J-1 

70  WRITE (6, 6)  J1,R(J) ,D(J) ,E(J) 

100  CONTINUE 
GOTO  10 
999  STOP 

1 FORMAT (212) 

2 FORMATCQ39.32) 

3 FORMAT Cl’, 'VALUES  OF  THE  POLYNOMIALS  IN  THE  RATIONAL  APPROXIMATIO 

:N  OF  1F1(AP;CP;-Z) '//'  ',’N-  ' .I2.T20, *AP  - '.Q39.32//'  ',T20,'CP 
: - '.Q39.32//'  Z-  ',Q39.32//'  I' ,T24, 'A(I) ' ,T65, 
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: 'B(I)  /) 

4 FORMAT ( ' ,I2,2X,Q39.32,2X,Q39.32) 

5 FORMAT (' O' VALUES  OF  THE  APPROXIMATION,  1ST  DIFFERENCES  AND  APPRO 
: XIMATE  ERRORS'//'  I',T12,'I-TH  APPROXIMATION  — F(I)',T47, 

: '1ST  DIFF' 'S. ' ,T60, ‘F(N)-F(I) '/) 

6 FORMAT ( ' ’ ,I2,2X,Q39.32,2X,Q10.3,2X,Q10.3) 

END 


C * THIS  SUBROUTINE  RETURNS  VALUES  A(I)  AND  B(I) , 1-1 N+l  * 

C * OF  THE  NUMERATOR  AND  DENOMINATOR  POLYNOMIALS  IN  THE  RATIONAL  * 

C * APPROXIMATION  OF  1F1(  AP  ; CP  ; -Z  ) . * 

C * * 

C * NO  OTHER  SUBROUTINES  ARE  CALLED  BY  THIS  ONE.  * 

C 1rk*1ckickickickirk-)ck-k*ick-kkk**-\kkkX+kkki'1t+k1ck*iticki'i'**i'*ir*irh**k1e+ie+kirk1e-kk*. 

IMPLICIT  REAL*I6 (A-H , O-Z) 

DIMENSION  A(l) , B(I) 

DATA  ZERO/O. QO/, ONE/1. QO/, TWO/2. QO/, THREE/3. QO/ 

C 

C INITIALIZATION  : 

C 

CT1-AP*Z/CP 
XN3-ZER0 
XN1-TWO 
Z2-Z/TWO 
CT2-Z2/ (ONE+CP) 

XN2-0NE 

A(l)-ONE 

B(l)-ONE 

B(2)-0NE+(0NE+AP)*Z2/CP 
A(2)-B(2) -CT1 

B ( 3) -ONE+(TWO+B ( 2) ) * (TWOfAP) /THREE*CT2 
A(3)-B(3) - ( ONE+CT 2 ) *CT 1 
CT1-THREE 
XNO" THREE 
C 

C FOR  1-3 N , THE  VALUES  A(I+1)  AND  B(I+1)  ARE  CALCULATED 

C USING  THE  RECURRENCE  RELATIONS  BELOW. 

C 

DO  100  1-3, N 
C 

C CALCULATION  OF  THE  MULTIPLIERS  FOR  THE  RECURSION 
C 

CT2-Z2/CT1/ (CP+XN1) 

Gl-ONE+CT2*(XN2-AP) 

CT2-CT2*(AP+XN1)  / (CP+XN2) 

G2-CT2*((CP=XN1)+(AP+XN0)  / (CT1+TW0)*Z2) 

G3-CT2*Z2*Z2/CT1/ (CT1-TW0)*(AP+XN2) / (CP+XN3)*(AP-XN2) 

C 

C THE  RECURRENCE  RELATIONS  FOR  A(I+1)  AND  B(I+1)  ARE  AS  FOLLOWS 
C - 

c 

B(I+l)“Gl*B(I)+G2*B(I-l)+G3*B(I-2) 

A(I+l)=Gl*A(I)+G2*A(I-l)+G3*A<I-2) 

C 

XN3-XN2 

XN2-XN1 

XN1-XN0 

xno-xno+one 

100  CT1-CT1+TV70 
RETURN 
END 


VALUES  OF  THE  POLYNOMIALS  IN  THE  PATIONAL  APPROXIMATION  OF  IFI(AP;CP; 
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THIS  MAINLINE  PROGRAM  ILLUSTRATES  THE  SYNTAX  OF  THE  SUB-  J 
ROUTINE  ' CIF1 ' WHEN  USED  TO  GENERATE  COEFFICIENTS  IN  THE  POLY-  .! 
NOMIALS  FOR  THE  RATIONAL  APPROXIMATION  OF  1F1 (AP ; CP ; -Z) . 


IMPLICIT  REAL*16(A-H,0-Z) 
DIMENSION  CA(26) tCB(26) ,NO(25) 
READ (5, 2, END-99 9)  AP.CP 
WRITE (6, 3)  AP.CP 
READ(5 , 1)  M, (NO(J) , J-1,M) 

DO  100  I-1,M 
N-NO(I) 


CALL  C1F1 (AP , CP , CA , CB , N) 


IN  THE  ABOVE: 

AP  IS  THE  NUMERATOR  PARAMETER  OF  THE  1F1 

CP  IS  THE  DENOMINATOR  PARAMETER  OF  THE  1F1 

N IS  THE  DEGREE  OF  THE  POLYNOMIALS  IN  THE  RATIONAL 

APPROXIMATION 

CA  AND  CB  WILL  CONTAIN  THE  COEFFICIENTS  IN  THE  NUMERATOR  AND 
DENOMINATOR  POLYNOMIALS  RESPECTIVELY 


NOTE  : THE  COEFFICIENTS  OF  THE  K-TH  POWER  OF  Z WILL  BE  PLACED 
IN  CA(K+1)  OR  CB(K+1)  AS  APPROPRIATE 


Nl-N+1 

100  WRITE(6,4)  N, (CA(J) ,CB(J) .J-l.Nl) 

GOTO  10 
999  STOP 

1 FORMAT (2612) 

2 FORMAT (Q39. 32) 

3 FORMAT (' 1' , 'COEFFICIENTS  FOR  THE  RATIONAL  APPROXIMATION  OF  1F1(  AP 
: ; CP  ; -Z  )'//'  \T20,'AP  - '.Q39.32/'  ',T20,'CP  - \Q39.32/) 

4 FORMAT ( ' ','N-  ' , I2.T18 , ' CA(I)  ' ,T58 , 'CB(I)  ' // 

: 26 (1X.Q39 . 32 , 2X, Q39 . 32/) /) 

END 


noon  o oooo  non 
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C 

C 

c 

c 

c 

c 

c 

c 


SUBROUTINE  C1F1(AP,CP,A,B,N) 


* THIS  SUBROUTINE  RETURNS  COEFFICIENTS  A(I)  AND  B(I)  , 

* I - 1,2, . . . ,N+1  , OF  THE  NUMERATOR  AND  DENOMINATOR  POLYNOMIALS 

* RESPECTIVELY,  IN  THE  RATIONAL  APPROXIMATION  OF  ORDER  N FOR 

* 1FI(  AP  ; CP  ; -Z  ). 

* 


* NO  OTHER  SUBROUTINES  ARE  CALLED  BY  THIS  ONE. 

**************  ******  *»*»★*■  A AAAAA*A**  *#**  Ah  A A A***  Hr***  ******  ********* 


IMPLICIT  REAL*16 (A-H , O-Z) 
DIMENSION  A(1),B(1) 

DATA  ONE/1. QO/, ZERO/O. QO/ 


INITIALIZATION  : 

XN-N 
XN1I-XN 
CP1-CP-ONE 
B(l) -ONE 
A(l)-ONE 
XI -ONE 
XIJ-ZERO 
DO  100  I-l.N 
Il-I+l 


FOR  I - 1,2 N , B(I+1)  IS  COMPUTED  AS  FOLLOWS 


B(I1)-(AP+XN1I) / (CP1+XN1I)*XN1I/ (XN+XN1I)*B(I) /XI 

A(I1)-0NE 
DO  50  J-1,I 

TO  CALCULATE  A(I+1) , WE  EMPLOY  B(J)  , J - l"2"”"”l+rAs"roLLOWS 


A(I1)-B(J+1) - (AP+XIJ) / (CP+XIJ)*A(I1) / (ONE+XIJ) 

50  XIJ-XIJ-ONE 

XIJ-XI 
XN1I-XN-XI 
100  XI-XI+ONE 
RETURN 
END 


****** 
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THIS  MAINLINE  PROGRAM  ILLUSTRATES  THE  SYNTAX  OF  THE  SUB- 
ROUTINE ' R1F1P*  FOP.  GENERATING  VALUES  OF  THE  NUMERATOR  AND 
DENOMINATOR  POLYNOMIALS  IN  THE  PADE  APPROXIMATION  OF 
1F1(  1 ; CP  ; -Z  ) . 


IMPLICIT  REAL*16 (A-H , O-Z) 

DIMENSION  A(26) ,B(26) ,R(26) ,D(26) ,E(26) 
DATA  ZERO/O.QO/ 

READ (5,1, END-999)  N , M 
READ (5, 2)  CP 
Nl-N+1 
D(N1)-ZER0 
E(Nl)«ZERO 
DO  100  I— 1 ,M 
READ(5 , 2)  Z 


CALL  RlFlP (CP , Z , A , B , N) 


IN  THE  ABOVE  : 

CP  IS  THE  DENOMINATOR  PARAMETER  OF  THE  1F1 

Z IS  THE  VALUE  OF  THE  ARGUMENT 

A AND  B WILL  CONTAIN  THE  VALUES  OF  THE  NUMERATOR  AND  DENOMINATOR 
POLYNOMIALS,  RESPECTIVELY,  FOR  ALL  DECREES  FROM  0 TO 
N INCLUSIVE 

N IS  THE  MAXIMUM  DEGREE  FOR  WHICH  VALUES  OF  THE  POLYNOMIALS 

ARE  TO  BE  CALCULATED 

NOTE  : VALUES  OF  THE  K-TH  DEGREE  POLYNOMIALS  WILL  BE  PLACED  IN 
A(K+1)  AND  B(K+1)  RESPECTIVELY. 


R(N1)-A(N1)/B(N1) 

DO  50  J-l.N 
Jl-Nl-J 

R(J1)«A(J1)/B(J1) 

D(J1)-R(J1+1)  -R(J1) 

50  E(J1)-R(N1) -R(J1) 

WRITE(6,3)  N , CP , Z 
DO  60  J-l.Nl 
Jl-J-1 

60  WRITE  (6 ,4)  J1,A(J),B(J) 

WRITE(6,5) 

DO  70  J-l.Nl 
Jl-J-1 

70  WRITE (6 ,6)  J1,R(J) ,D(J) ,E(J) 

100  CONTINUE 
GOTO  10 
999  STOP 

1 FORMAT (212) 

2 FORMAT (Q39. 32) 

3 FORMAT ('l',' VALUES  OF  THE  POLYNOMIALS  IN  THE  PADE  APPROXIMATION  OF 

: 1F1 (1 ; CP ; -Z) * / / * ' , 'N  - ' , 12 ,T20 , ’CP  - '.Q39.32 //’  Z-  1 

:Q39. 32// ' I' ,T24, 'A(I) ' ,T65, 'B(I) ' /) 

4 FORMAT ( ' ' ,I2,2X,Q39.32,2X,Q39.32) 
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rOR!«AT('0'  , ' VALUES  OF  THE  APPROXIMATION,  1ST  DIFFERENCES  AND  APPRO 
:XD1ATE  ERRORS  * / / * I',T12,'I-TH  APPPOXIFATION  — F(I) ' , 

:T47 , ' 1ST  DIFF"S . \T60 , 'F(N) -F(I)  ' /) 

FOR! ’AT ( ' ' ,I2,2X,Q39.32,2X,Q10.3,2X,Q10.3) 

END 


noon  non  o o o o non 
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C 

C 

C 

C 

C 

C 

C 


* THIS  SUBROUTINE  RETURNS  VALUES  A (I)  AND  B(I) , I-1,...,W-1  * 

* OF  THE  NUI  IEP.ATOR  AND  DENO!  TIN  ATOP  POLYNOMIALS  IN  THE  PADE  * 

* APPROXIMATION  OF  1FI(  1 ; CP  ; -Z  ) . * 


IMPLICIT  REAL*16  (A-11 , 0-2) 
DIMENSION  A(I)  ,B(1) 

DATA  ONE/1.QO/,TUO/2.QO/ 


INITIALIZATION  : 


XI1-0NE 

B(1)-0KE 

A(1)-0NE 

CT1-CP+0NE 

CP1-CP-01IE 

ZZmZ*Z 

B ( 2)  -ONL+Z  / CT1 
A(2)-B(2)-Z/CP 

FOR  1-2 N , THE  VALUES  A<I+1)  AND  B(I+I)  ARE  CALCULATED 

USING  TIIE  RECURRENCE  RELATIONS  BELOW. 

DO  100  1-2, N 

CALCULATION  OF  THE  MULTIPLIERS  FOR  THE  RECURSION 
CT2«CT1*CT1 

G1-0NE+CP1/ (CT2+CT1+CT1)*Z 
C2-XI1/  (CT2-0NE)*(XI1+CP1) /CT2*ZZ 


THE  RECURRENCE  RELATIONS  FOR  A(I+1)  AND  B(I+1)  ARE  AS  FOLLOWS 


A(I+1)-G1*A(I)+G2*A(I-1) 
B(I+1)-C1*B(I)+G2*B(I-1) 
C 

CT1-CT1+TU0 
100  XI 1— XI  l+ONE 

RLTURN 
END 


VALUES  OF  THE  POLYNOMIALS  IN  T1IE  PADE  APPROXIMATION  OF  1F1(1;CP;-Z) 
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THIS  MAINLINE  PROGRAM  ILLUSTRATES  THE  SYNTAX  OF  T1IE  SUB- 
ROUTINE 1 C1F1P ' WHEN  USED  TO  GENERATE  COEFFICIENTS  IN  THE 
POLYNOMIALS  FOR  THE  PADE  APPROXIMATION  OF  1F1(1;CP;-Z) . 


IITPLICIT  REAL*16  (A-H , O-Z) 
DECISION  CA(26),CB(26),NO(25) 
READ(5, 2, END-999)  CP 
WRITE (6, 3)  CP 
RLAD(5,1)  K.(NO(J),J-l,M) 

DO  100  I-l.M 
N-NO(I) 


CALL  C1F1P (CP , CA , CB , II) 


IN  THE  ABOVE: 


CP  IS  THE  DENOMINATOR  PARAMETER  OF  THE  1F1  IN  THE  PADE 

APPROXIMATION 

N IS  THE  DEGREE  OF  THE  POLYNOMIALS  IN  THE  PADE 

APPROXIMATION 

CA  AND  CB  WILL  CONTAIN  THE  COEFFICIENTS  IN  THE  NUMERATOR  AND 
DENOMINATOR  POLYNOMIALS  , RESPECTIVELY 


NOTE  : THE  COEFFICIENTS  OF  THE  K-TH  POUEP  OF  Z WILL  BE  PLACED 
IN  CA(K+1)  OR  CB(K+1)  AS  APPROPRIATE 


100 

999 

1 

2 

3 


Nl-N+1 

I HITE (6, 4)  N, (CA(J) ,CB(J) , J-1.H1) 

GOTO  10 
STOP 

FORMAT (2612) 

FORI  AT  (Q  39 . 3 2 ) 

FORMAT ( ' 1 ' ' COEFFICIENTS  FOP  THE  PADE  APPROXIMATION  OF  IF1(  1 ; 
: ; -2  )’//‘  \T20,'CP  - ' ,039. 32/) 

FORMAT ('  \'N  - ' ,I2,TI8,  CA(I) ' ,T50, 'CB(I) '// 
:26(1X,039.32.2X,Q39.32/)/) 

END 


CP 


nano  n nnnn  non 
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SUBROUTINE  C1F1P(CP,A,B,N) 

C ***■**'*****-**■**■*'*****  A A*  **************  Irltitltieicitlrklcivlt-kiFkiilrk  AnAArAArA*  ***** 

C * THIS  SUBROUTINE  RETURNS  COEFFICIENTS  A(I)  AND  B(I)  , 

C * I - 1.2 N+l  , OF  THE  NUMERATOR  AND  DENOMINATOR  POLYNOMIALS 

C * RESPECTIVELY,  IN  THE  PADE  APPROXIMATION  OF  ORDER  N FOR 

C * 1F1 ( 1 ; CP  ; -Z  ). 

C * 

C * NO  OTHER  SUBROUTINES  ARE  CALLED  BY  THIS  ONE. 

C *•***■*■*■**"*  ******  A A A AAAArAAAArAA^*********1*  AAA  A*-*****  '*"***  At**  A A*  A AAA  A AAA 

IMPLICIT  REAL*16(A-H,0-Z) 

DEIENSIOII  A(l)  ,B(1) 

DATA  ONE/1. QO/, ZERO/O. QO/ 

INITIALIZATION  : 

XN-N 
XN1I-XN 
B(l)"ONE 
A(l)-ONE 
XI-ONE 

CP2NI-CP+XN+XN-ONE 
XIJ-ZERO 
DO  100  I-1,N 
Il-I+l 


FOR  I - 1,2, ...  ,H  , B(I-fl)  IS  COMPUTED  AS  FOLLOWS 


B(I1)-XN1I/CP2NI*B(I)/XI 

A(I1)-0NE 
DO  50  J-1,1 

TO  CALCULATE  A(I+1)  , WE  EMPLOY  B(J) , J - 1,2,... ,1+1  AS  FOLLOWS 


50 

100 


A(I1)-B(J+1)-A(I1)/(CP+XIJ) 

XIJ-XIJ-ONE 

XIJ-XI 

XN1I-XN-XI 

CP2NI-CP2NI-ONE 

XI-XI+ONE 

RETURN 

END 


* * * * * * 


COEFFICIENTS  FOR  THE  FADE  APPROXIMATION  OF  1F1(  1 
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